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Aims. During the first 40 s after their birth, proto-neutron stars are expected to be subject to at least two types of instability. The first 

(~| , one, the convective instability, is excited in the inner regions, where the entropy gradient produces a Rayleigh-type convection. The 

^ |. second one, the neutron-finger instability, is instead excited in the outer layers where the lepton gradients are large. Both instabilities 

I ' involve convective motions and hence can trigger dynamo actions that may be responsible for the large magnetic fields in neutron 

O ' stars and magnetars. However, because they have rather different mean turbulent velocities, they are also likely to give rise to different 

types of dynamo. 

^2 ■ Methods. We have solved the mean-field induction equation in a simplified one-dimensional model of both the convective and the 

^ ' neutron-finger instability zones. Although very idealized, the model includes the nonlinearities introduced by the feedback processes 

that tend to saturate the growth of the magnetic field (a-quenching) and suppress its turbulent diffusion (77-quenching). The possibility 

p^ \ of a dynamo action is studied within a dynamical model of turbulent diffusivity where the boundary of the unstable zone is allowed to 

►^ . move. A large number of numerical simulations have been performed in which the relevant parameters, such as the spin-period, the 

C^ ' strength of the differential rotation, the intensity of the initial magnetic field, and the extent of the neutron finger instability zone, have 

^j< ' been suitably varied. 

— ^ , Results. We show that the dynamo action can also be operative within a dynamical model of turbulent diffusivity and that the 

amplification of the magnetic field can still be very effective. Furthermore, we confirm the existence of a critical spin-period, below 
which the dynamo is always excited independently of the degree of differential rotation, and whose value is related to the size of the 
neutron-finger instability zone. We provide a relation for the intensity of the final field as a function of the spin of the star and of its 
differential rotation. 

C ■ Conclusions. Although they were obtained by using a toy model, we expect that our results are able to capture the qualitative and 

^^ ■ asymptotic behaviour of a mean-field dynamo action developing in the neutron-finger instability zone. Overall, we find that such a 

dynamo is very efficient in producing magnetic fields well above equipartition, and thus that it could represent a possible explanation 
for the large surface magnetic fields observed in neutron stars. 
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1 . Introduction rate is Bd ~ 2 - 4 x 10'^ G (Pavlov et aL 120021 ). while the sur- 

^, , ,• ^ , , , , face field estimated from the absorption features in its spectiTim 

The present understanding of the processes that produce the -^ g^^^ ^ j 5 ^ ^qu q (g^^^^j ^^ ^j 2002 ). Similarly, observa- 

large magiietic field strengths observed in neutron stars (NSs) ^.^^^ ^^ ^^^ j^^^ ^^^ B1821-24 (Becker et al.|2003) indicate 

IS still far from being complete. Most of the information about j. j. □ ^ in^ r while R ~ 10" T 
their magnetic fields is in fact derived either from their X-ray 

spectra, or from their spin-down when these NSs are seen as The existence of magnetic fields with diff^erent strengths and 
pulsars. While the former reflects a measure of the local surface distributed on different length-scales can be explained in terms 
field B.urf, the latter provides information on the global dipo- of a dynamo mechanism driven by the simultaneous presence 
lar magnetic field Bj, if the spin-down is assumed to be solely of rotation and turbulent motions. During the first ~ 40 s after 
due to dipolar electromagnetic emission. These two measures their birth, proto-neutron stars (PNSs) are expected to develop 
are not always in agreement, showing that the measured mag- hydrodynamical inst abiliti es (Epstein 1979; Livio et al. 19801 
netic fields may have very different length-scales and intensities. Burrows & Lattimer [198611, which can excite a hydro-magnetic 
In particular, they seem to suggest the presence of more intense dynamo. Such instabilities could be essentially of two types. The 
and smaU-scale surface magnetic fields, together with less strong first one, driven by the entropy gradient, is a Rayleigh-type con- 
globally dipolar ones. For the pulsar IE 1207.4-5209, for in- vective instability (CI) that operates in the inner regions of the 
stance, the dipolar magnetic field estimated from the spin-down star The second one is a double diff^usive instability, driven by 

both the entropy and leptonic gradients. This is usually referred 

Send offprint requests to: L. Naso to as the neutron-finger instability (NFI); it operates in the outer 
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regions of the PNS (Miralles et al. |2000t |2002] l and is expected 
to evolve by creating finger-like downflows when the neutrinos 
are still confined (Mezzacappa et al. 1998). Although some au- 
thors have recently raised doubts about the existence of the NFI 
(Buras et al. l2006l Dessart et al. ,2006. Bruenn et al. |2007] |. no 
firm conclusion has yet been reached, leaving the debate open. 
Here, we do not attempt to enter this debate but rather, because 
of the interesting astrophysical implications that it may have, we 
will consider the NFI as taking place and having the dynamical 
properties as described by Miralles et al. ( 120001 |2002] l. 

We note that the co-existence of the two instability mech- 
anisms produces both a local dynamo process (Thompson & 
Duncan[l993 , Xu & Busse 2001 ) and a mean-field one (Bonanno 
et al. I2003J . As shown in Miralles et al. (|2000l |2U02] i. the 
growth-times of the instabilities in the two regions differ by 2 
or 3 orders of magnitudes, being tqi ~ 0.1 ms in the CI zone 
and tnpj ~ 30 - 100 ms in the NFI zone. Since the typical spin 
period of a PNS is P ~ 100 ms, the turbulent eddies created by 
the CI are not influenced by the rotation and therefore they can 
only excite a local dynamo. On the other hand, the Rossby num- 
ber in the NFI zone, defined as Rq - P/tnfi, is about unity and 
the turbulent motions can therefore be influenced significantly 
by the rotation, favoring the excitation of a global mean-field 
dynamo. Because of the large difference in the growth-rates of 
the two instabilities, the two processes, i.e. the local dynamo and 
the global mean-field one, are essentially decoupled. 

Here, we focus our attention on the turbulent mean-field dy- 
namo action that may be excited in the NFI zone. More specifi- 
cally, we exploit a simple one-dimensional (ID) toy model that 
aims at capturing, at least qualitatively, the features of the dy- 
namo action. The model, which includes the nonlinearities in- 
troduced by the feedback processes, which in turn tend to satu- 
rate the growth of the magnetic field, i.e. a-quenching (Bonanno 
et aI. l20U5ll2006l Rudiger & Arh rT996] l, and suppress its turbu- 
lent diffusion, i.e. ?7-quenching (Rudiger & Arlt 1996), is evolved 
numerically with a very large variety of initial conditions. These 
include varying the spin period of the PNS, the strength of the 
differential rotation between the core and the surface, the inten- 
sity of the primordial (seed) magnetic field, and the extent of the 
NFI zone. 

Overall, we find that increasing the extent of the neutron- 
finger instability zone favours the dynamo excitation, and that 
the combined action of differential rotation and diffusion can 
produce an increase in the strength of the generated toroidal 
field by several orders of magnitude. We also confirm the ex- 
istence of a critical spin-period, below which the dynamo is al- 
ways excited independently of the differential rotation strength, 
and whose value is related only to the size of the neutron-finger 
instability zone. Since the numerical simulations show that the 
turbulent mean-field dynamo is very efficient in producing mag- 
netic fields well above the equipartition value, we confirm the 
results of Bonanno et al. '2006' even in presence of both a mov- 
ing boundary of the instability zones and of an ?7-quenching, i.e. 
the proposed mechanism could represent an explanation for the 
large surface magnetic fields observed in several neutron stars. 

The paper is organized as follows: in Sect. |2]we formulate 
the model and derive the system of partial differential equations, 
which are then solved with the numerical approach described 
in Sect. [3] In Sect. |4] we describe the methodology and report 
numerical results. Finally in the conclusions we summarize our 
results. 




Fig. 1. Schematic diagram of our model of the PNS. 



2. The model 

The mean-field induction equation of a turbulent, magnetized 
and conducting plasma can be written as 



5rB = V X (v X B -H aB) - V X (77V X B) 



(1) 



where B is the mean magnetic field, v the mean velocity field, t] 
the turbulent magnetic diffusivity, and a a pseudo-scalar measur- 
ing the efficiency of the dynamo a-effect (this is usually referred 
to as the a-parameter). This effect is related to the generation of 
a net electromotive force &nirb by the turbulent components of 
the magnetic field (B') and the velocity field (v')- More specif- 
ically, &„irh = v' X B' = CKijBj + PijitdkSj -H ... In the case of 
isotropic turbulence a/j = aSij and /Sjjk = /Se/jk and S,uf-h can be 
approximated as &nirb = aB - j3V x B. 

Our simplified model follows closely the one introduced by 
Brandenburg et al. (119891) . Rudiger et al. (1994), and Rudiger 
& Arlt ( 1996). It uses orthogonal Cartesian coordinates, and the 
PNS is modeled has having planar symmetry with the z-axis be- 
ing the axis of rotation, the (x, 3')-plane representing the stellar 
equatorial plane, and H the semi-height being the radius of the 
PNS. 

Clearly, this simplified model has the advantage of leading to 
a very simple expression for the induction equation ([T]l which, 
after making a suitable choice for the velocity and magnetic 
fields, can be recast into a system of two coupled partial differ- 
ential equations that are of first order in time. This reduces the 
computational costs enormously and allows a parametric investi- 
gation to be performed, which would be impossible if it involved 
fully 3D simulations. Yet, despite the considerable simplifica- 
tions, this simple geometry has been shown to yield instructive 
results that have opened the way for modern, realistic 3D models 
(Weiss et al. 119841 ). Our simplified model of the PNS is shown 
in Fig.[ri where we highlight the CI and NFI zones that are sep- 
arated by a thin interface of thickness d, which is not shown in 
the figure. Note that despite the cylindrical appearance of the 
star, the system does not have a cylindrical symmetry but rather 
a planar one across the (x,y) plane. 

Another important aspect in the present modeling of the tur- 
bulent mean-field dynamo is the so-called kinematic dynamo 
approximation, in which the solution of the induction Eq. ([T]) 
is assumed to be decoupled from the Euler equations (i.e. the 
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PNS is assumed to be in hydrostatic equilibrium at all times), 
and no feedback from the magnetic field is taken into account in 
the conservation of momentum and energy. As a result, the ve- 
locity field is taken as pre-assigned and time-independent. This 
approximation works quite well as long as the magnetic field 
strength is small (i.e. as long as the magnetic pressure and ten- 
sion are negligible in the Euler equations), but has the drawback 
that nothing prevents the indefinite growth of the magnetic field 
once a dynamo action is present. In real systems, the velocity 
profile will adjust itself in such a way as to reduce the efficiency 
of the dynamo and here, as a way of mimiking this feedback, we 
introduce two "quenching functions" (described later in more 
detail) that suppress the amplification as the total field increases 
above a certain threshold, thus leading the system towards satu- 
ration even in the absence of a consistent feedback. 

Our model for the kinematic dynamo follows the one pro- 
posed by Riidiger et al. ( |19941 l and Blackman & Brandenburg 
(I2OO2I 1 and, in particular, we consider a velocity field 

V = (0, kx, 0) , 

and a magnetic field having components only in the (x, y) plane 
but a with a vertical dependence 



B = {bAz),B,.(z),0) 



The a-parameter and the magnetic diffusivity t] are both ex- 
pressed as the product of three terms 



a - ao a'(z,t) if/J^) , 
ri ^ m r]'(z,t)i//,,{B) , 



(2) 
(3) 



where ao and 770 measure the strength of the a-effect and turbu- 
lent diffusion respectively, while a'(z, t) and rj'iz, t) represent the 
profiles of a and rf in the two instability regions. More specif- 
ically, a' is chosen to be antisymmetric across the equatorial 
plane and different from zero only in the NFI zone, where the 
mean-field dynamo is at work. The turbulent diffusivity rj' , on 
the other hand, is set to be of the order of unity in the NFI zone 
and about an order of magnitude larger in the turbulent CI zone. 
We implement these prescriptions by making use of the error 
function "erf" as 



a\z, t) = 



i[l+erf(-(z + i)/flO] 



[1-Herf((z-^)/^] 



z^{-H,0] 



z&[0,H] 



(4) 



i{z,f)^-{\0- 



l+&rf[{z-X)ld] 



10-- l+trf[-{z + X)ld] 



(5) 



where A= z\ + Vt, with z\ being the coordinate of the boundary 
between the CI and NFI zones and V the expansion velocity of 
the boundary layer The quantity d represents the thickness of the 
interface between the CI and NFI zones and is used to obtain a 
smooth change of the error function, with smaller values leading 
to sharper changes; for the results reported here we have chosen 
d/H = 0.04. The profiles of a' and 77', as given by Eqs. (|4|i and 
(|5]), are shown in Fig. |2] with the solid and dashed lines indi- 
cating the initial conditions and the final conditions after 40 s, 
respectively. 

Finally, ^^^(B) and t//^{R) appearing in Eqs. (|2]) and ([3]) rep- 
resent the quenching functions for the c-effect and turbulent dif- 
fusion rj, respectively. These terms are used to limit the oth- 
erwise unlimited growth of the magnetic field (a-quenching) 



1 




^C ' 


: \ \ : 


-1 


\ \ : 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-0.5 0.5 1 



10 


1 




•■ 1 \ •■ - 


^'5 


- 


'-- 




_ 
/ 


I ■ 


n 





-0.5 



0.5 1 



Fig. 2. Top panel: profiles of the normalized function a' at the 
beginning of the numerical simulation for f = (solid line) and 
at the end of the numerical simulation for f = 40 s (dashed line); 
bottom panel; the same as in the top panel but for 77'. 



and suppress its turbulent diffusion (//-quenching). In general 
they are expected to have a different dependence on the mag- 
netic field strength, but we here consider a single expression for 
the two functions in terms of the equipartition magnetic field 
Bpn - (u) -\J4-Kp, where (u) is the mean velocity of turbulent ed- 



dies and p the mass-density 



tA«,,(B) = 



1 + 



J-H\Beq 



dz 



(6) 



Making use of these definitions, Eq. ([T]i can be split into the 
following two coupled scalar partial differential equations 



drB, = -d,(aB,) + d, (T]d,B,) , 
d,By = d,(aB,) + B,d,Vy + d- [rjdrB,) 



(7) 
(8) 



which can also be written in a dimensionless form by scaling 
lengths in units of the semi-height of the cylinder H, times in 
units of the diffusion time t^ - H^/tjq, and magnetic fields in 
units of Beq- We also find it useful to introduce the dimensionless 
parameters 

r - H 

Ca = ao— , 



and Cn = <9^v,. — \^k\ — 

/no) \m 



(9) 



which represent the Reynolds numbers for the a-effect and the 
differential rotation, respectively. Furthermore, by introducing 
the standard vector potential A = [Ax(z),Ay{z),Q], so that the 
poloidal component of the magnetic field is Sp = B^ = -d^Ay, 
we obtain the following dimensionless equations for the vector 
potential J?l - Ay and the toroidal component of the magnetic 
field S, = B, 



d,J[ = CMz, t)ilf,(B,o,)S, + i](z, t)ilf,(B,o,)dt/^ 



(10) 
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d,S, = -Cad, [a(z, t)iff„(B,o,)d,m] - Cnd.Jl + 
+ d,[j](z,t)ilr^(B,o,)d,S,] 



(11) 



where B,o, = [S^+(5:J?l)^]'^l Once the initial conditions S,(z,0) 
and y[(z, 0) are given, together with the parameters, €„, Cn, a, 
Tj, iffa, if/,j, and suitable boundary conditions at the stellar edges, 
it is possible to solve Eqs. ( fTOb and ( fTTI ) to describe the time evo- 
lution of the magnetic field. Our choice for the boundary condi- 
tions reflects the fact that we are interested in an adiabatic evo- 
lution of the magnetic field in the stellar interior, thus neglecting 
the energy losses related to a Poynting flux. Because of this, we 
simply set the magnetic fields at these locations to zero. 

It is important to note that the parameters €„ and Cn are not 
linearly independent but can be related in terms of the strength 
of the differential rotation, of the pressure scale-height and of 
the "radius" of the PNS. To deduce this relation we recall that 
the aQ-dynamo assumes that ao ~ ilLp, where Q is the angular 
velocity and Lp the pressure scale-height, so that 



P 27tL„H 
P ~ -^ ^ P— 

Ca CaTJO 



(12) 



where ^ = InLpH/rjo. For a typical PNS with mass ~ IM©, 
H ~ 15 km and Lp ~ 3 km (Bonanno et al. l2003l l. t/o ~ Lp/T"NFI 
ranges from 9 x 10" cm^s"' for NFI eddy turnover times of 
100 ms, up to 3 X lO'^cm^s"' for NFI eddy turnover times of 
30 ms. Such values of tjq yield typical diffusion timescales r^ 
ranging from 0.75 s to 2.5 s, respectively. As a result, the param- 
eter ^ is expected to be roughly in the range 1 s < ^ < 3 s for all 
of the relevant parameter space considered here. 

In a similar way, since 5 vV,, = A: ~ AQ, and defining the rela- 
tive differential rotation strength as ^ = AO/O ^ (Qj - i^c)/fis, 
where Q, and Qc are the angular velocities of the surface and the 
core respectively, it is possible to conclude that 






(13) 



10'^ gem ^ and that the eddy convective velocities (u) ^ 
it follows that Bgq is of the order of 



where ( = Lp/H ^1/5. Here we use AQ as a global measure of 
the rotational stress instead of using a necessarily arbitrary func- 
tion that describes the behaviour of Q(z) in the region between 
fie and Qj. On the other hand we are severely limited by our ig- 
norance of the detailed processes leading to the appearence of 
the differential rotation in the NFI zone of PNSs. 

For all of the calculations reported here we assume that the 
angular velocity of the core is larger than that of the surface, so 
that ^ < and, conservatively, we limit our analysis to values 
of 1^1 not exceeding 10^, i.e. a core rotating 10^ times faster than 
surface. Finally, by assuming a mass-density in the NFI zone 

P 

Lp/TNFI ~ 3 X 10^ cms 

lO'^ G. The parameters defined in Eqs. (fTST i and ( fTST l essentially 

determine the parameter space for the solutions of Eqs. ( fTOl l and 

(ED. 

3. Numerical method and tests 

In order to solve the mixed parabolic -hyperbolic system of par- 
tial differential equations ( fTOl i and ( fTTT i. we discretize the con- 
tinuum space-time by replacing it with a two dimensional grid, 
where the two dimensions represent the space and the time vari- 
ables, z and t, respectively. We use constant spacing in both di- 
rections, with a typical grid of 50 zones. Tests were performed 
with a larger number of gridpoints (100, 200 and 400) and have 



revealed that a minimum of 50 gridpoints was sufficient to yield 
a small-enough truncation error The evolution algorithm cho- 
sen is the FTCS (Forward-in-Time, Centered-in-Space) scheme, 
which gives a first-order approximation for the time deriva- 
tives and a second-order approximation for the space derivatives. 
Furthermore, stability requires the timestep to be At = O(Az^) 
(we typically use At = 10"^ Az-), thus making the whole algo- 
rithm second-order both in space and in time. Using ghost-zones 
for implementing the boundary conditions, the final form of the 
finite-difference equations is 



Jl" 



^"j+Atd"jS"j 



on+l 
J 



where 



yj-- 



and 






rzin (Tin 

2Az 



(14) 



Az2 



k?^ 



Caa"jK 



At 
2Ai 
b^-Cn 



d" 
J 



C — i^Qr 



2Az 



(15) 

(16) 

(17) 
(18) 

(19) 



As a test of the code we have considered the equations when 
Ca = Cq = 0, ij/^ - I and rjiz, f) = 5 ■ 10"^ so as to have two 
decoupled purely parabolic equations, and compared the numer- 
ical solution with the analytic one. The result of the comparison 
is that the maximum error is of the order of 2%. 

We have also checked the convergence of the method by 
comparing the numerical solutions obtained with different spa- 
tial grids at a given time. The same continuum function u(z, t) 
is approximated in a different way according to the space-time 
discretization: u(z, t) - u"{h) + e/,, where e/, is the truncation 
error. If we suppose that this error depends on the chosen space- 
interval h only, we can write: e/, - k h'\ with p being the order of 
convergence. It is then possible to deduce the following relation 
among the values of different discretized functions calculated at 
the same grid location /: 



m;'(/i) - uiih/i) 

ulih/l) - u'!(h/4) 



2P 



(20) 



The value of p obtained as averaged over the spatial domain is 
2.00 + 0.01, thus demonstrating a second-order convergence. 



4. Analysis and results 

While this analysis aims at a better understanding of the be- 
haviour of the dynamos operating in the first stages of the life 
of a PNS, it is a long way from reproducing realistic condi- 
tions. This is partly due to the simplicity of the model employed, 
and partly to the still poorly constrained physical conditions of 
a newly born PNS. We recall that, according to Miralles et al. 
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Table 1. The parameters z\ (initial position of the boundary 
layer), V (mean velocity of the boundary layer), and 77-q (77- 
quenching activated or not) that define the configurations ana- 
lyzed, and the ranges of values of Bs, Ca and \q\ used for the 
simulations. As regards Ca, vve mainly used values between 5 
and 200, since for smaller values the dynamo is not excited and 
for larger values the spin period would be too short. 





z'l 


V 


'7-q 


BJB,^ 


Ca 


\q\ 


A 


0.3 


0.0 


no 


10-' - 10-' 


< w 


10-» - 10^ 


B 


0.8 


0.0 


no 


10-^ - 10-' 


< 10^ 


lO-** - 10^ 


Aq 


0.3 


0.0 


yes 


10-7 _ 10-1 


< 10^ 


lO-'* - 102 


Bq 


0.8 


0.0 


yes 


10-7 _ 10-1 


< 10^ 


lO-'* - 102 


AB 


0.3 


0.0125 


no 


10-' 


< 10" 


10-'' - 10" 


ABq 


0.3 


0.0125 


yes 


10-7 


< 10" 


lO-** - 10" 



([2OOOI |2002] |. the NFI is expected to last only about 40 s, dur- 
ing which the NFI zone goes from occupying a large fraction of 
the envelope, to being confined to a small layer and then disap- 
pearing completely. We model this by assuming that the initial 
position of the NFI-CI boundary layer is at z\ - 0.3// and the 

final one, after 40 s, is at z\ — 0.8// (cf. Fig.[lJ, with an average 

expansion velocity of the layer that is V = {z\ -z\)/t = 188 m s"' 
for // = 15 km. In this case we also find it convenient to express 
all variables in terms of dimenionless quantities after introducing 
z' = z/H, t' = t/T^, and V = Vt^/H. As a result, the coordi- 
nate position of the CI-NFI layer appearing in Eqs. (|4| and (|5]l 
can be written as /I = z\ + V'f - 0.3 -1- 0.0125 (s"') f (s), where 
-1 < z' < 1, and < ? < 40 s, corresponding to a number of 
diffusion times ranging from 16 to 53, depending on the turnover 
time of the NFI eddies. 



4.1. Initial models 

As a representative sample of initial data we have considered six 
different models, four of which have the size of the NFI zone 
being constant in time, with either z[ = 0.3 (large instability 
zone) or z[ - 0.8 (small instability zone), and which are there- 
fore referred to as static. For each of the two values of z'j we have 
examined the behaviour with and without 77-quenching. Besides 
the static configurations, which are useful for studying the ex- 
treme cases of thick and thin NFI zones respectively, we have 
also considered two cases, which are referred to as dynamical, 
in which the NFI is allowed to shrink in time from the initial 
value of z'j - 0.3 to the final one of z'j = 0.8, over the 40 s during 
which the instability is expected to be active. For these dynam- 
ical models we have also studied the effect of activating or not 
activating the 77-quenching. Such cases are useful for studying 
the role played by a dynamically shrinking NFI zone in the onset 
of the dynamo action and in the final magnetic fields obtained. 

The static models are indicated as A, B, Aq, Bq, while the 
dynamical ones are indicated as AB and ABq; in all cases, the 
letter "q" is used to indicate whether or not 77-quenching is taken 
into account. For all of these models we have carried out a large 
number of simulations by varying the seed magnetic field B^, 
as well as the spin period P and the strength of the differential 
rotation \q\; all the models have been evolved for 40 diff'usion 
times (i.e. between 30 s and 100 s). A summary of the properties 
of the different initial models and of the parameters used in the 
simulations is presented in Table [1] 



4.2. Time evolution and critical period 

The time evolutions of the average toroidal field S, and poloidal 
field Sp (calculated in terms of their 2-norms) are shown in Figs. 
|3] and Irrespectively for the configuration Aq, with \q\ - 2 and 
C„ = 4 (0.25 s < P < 0.75 s), and for diff'erent values of the seed 
magnetic field B^ in the range lO-^Bgq ^ Bs < lO-'figq. 

For all of the configurations examined, the evolution of the 
magnetic field is rather similar and can be separated into two 
main stages: the first one is a transient phase during which the 
dynamo action amplifies the seed magnetic field exponentially; 
in the second phase the magnetic field instead reaches saturation 
around the equipartition value through the back reaction of the 
ff-quenching. In general, the field reaches saturation within the 
40 s lifetime of the NFI zone, except when the seed magnetic 
field is lower than 10-^Beq and \q\ < 2. Nevertheless, in these 
cases the final magnetic value is of the order of 10-^Beq, thus 
corresponding to 10" G. 

These two stages can easily be distinguished in Figs. [3] and 
m which report the time-evolution of the toroidal and poloidal 
magnetic fields, respectively. The secular slope of the curves in 
the exponential-amplification phase is clearly independent of B^ 
and constant in time, with B - B^ e '^'^""'i' and Tamp ~ 0.6 td, 
where td - 25tnfi is the diff'usion timescale. The growth-time 
however depends on \q\ and Ca, thus suggesting that the seed 
magnetic field determines only the time interval necessary to 
achieve the saturation, but not the final strength of the magnetic 
field. Clearly, if the initial field is too small, the dynamo can- 
not reach the saturation phase within 40 s. Other values of the 
growth times are Ta,„p ~ 0.98 tq for configuration ABq (with 
C(j = 4 and \q\ - 3) and T„„p ~ 2.19 tq for configuration B 
(with Ca = 29 and \q\ = 4). 

It is possible that the combination of the rotation rate and 
differential rotation is not adequate to excite the dynamo and that 
the magnetic field is not amplified but rather decays with time. 
For each value of \q\ it is possible to find a threshold value of Co 
such that for higher rotation rates (shorter periods) the dynamo 
is excited, while for lower rotation rates (longer periods) it is not. 
This value of C„ defines a critical period Pc = Pdlql) through ^ 
[seeEq. (fntl. 

The presence of a critical period within the space of param- 
eters is shown in Fig.|5]for the configuration A. The behaviour 
of the dynamo solutions is plotted as a function of the PNS spin 
period P and diff'erential rotation strength \q\ for three different 
values of ^ (corresponding to different turnover times of the con- 
vective eddies in the NFI zone: ^ - I, tnfi - 30ms; ^ - 2, 
tnfi - 65 ms and f = 3, tnfi - 100 ms). The three curves repre- 
sent the thresholds between the regions in which the solutions of 
the dynamo equations grow in time (regions below the curves), 
thus allowing dynamo excitation, and the regions in which the 
solutions decay, thus preventing dynamo excitation. The three 
lines in Fig. |5] define therefore the critical dynamo period as a 
function of \q\, and it is evident that the critical spin velocity 
above which the dynamo operates decreases with increasing dif- 
ferential rotation. The behaviour of these curves suggests that is 
possible to define a global critical period Pc as the minimum of 
Pc, so that if a PNS is rotating with a period shorter than Pc then 
for that PNS the dynamo will be excited. 

The values of the global critical period obtained for the four 
static configurations are summarized in Table |2] Note that for 
a fixed value of ^ these periods depend only on the thickness 
of the NFI zone and are thus independent of the presence of 
77-quenching. As a result, in PNSs with larger NFI zones, the 
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Fig. 3. Behaviour of the 2-norm of "B, as a function of the num- 
ber of diffusion times {tItd) for the Aq configuration, \q\ - 2, 
Cq, = 4 (250 ms < P < 750 ms), and for different seed magnetic 
field strengths, Bs- The growth-rate of the magnetic field during 
its amplification phase, before reaching the saturation, is almost 
independent of the seed magnetic field. 



Table 2. For each static configuration the critical value of Ca 
(C^) and the corresponding global critical period of the PNS, 
Pc - ^/C^ ms, for ^ e [ 1 , 3] are reported. 
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Fig. 4. The same as in the Fig. |3]but for S^. 



dynamo action will be excited more easily, and thus at lower ro- 
tation rates, than for PNSs having smaller NFI zones. 



Fig. 5. Behaviour of the dynamo solutions as a function of the 
PNS spin period P and of the differential rotation strength \q\ 
for three different values of ^, related to the turnover time of the 
convective eddies in the NFI zone. 



4.3. The role of thickness of the NFI zone 

Since the NFI zone is the region where the mean field dynamo 
is at work, it is natural to expect that by increasing the size of 
this part of the star the efficiency of the dynamo will also in- 
crease. To quantify this "improved efficiency" we can compare 
the results obtained for two configurations: one in which the NFI 
zone occupies 70% of the star (configuration A) and another one 
in which it covers 20% (configuration B). Reducing the extent 
of the active part of the star has two main effects: reducing the 
critical period and reducing the intensity of the final magnetic 
field. 

The first of these effects has already been discussed in the 
previous section, and from Table |2] it is possible to note that 
reducing the NFI zone by a factor of 3.5 decreases the critical 
period by a factor of 6. To quantify the second effect we con- 
sider the ratio between the final strength of the magnetic field in 
configurations A and B, as a function of the differential rotation 
parameter \q\ and of Ca- As expected, the toroidal magnetic field 
for the configuration A is larger than the one for configuration 
B, regardless of any other parameter The value of the ratio is re- 
ported in Table [3] and it should be noted that it is really only the 
toroidal-field ratio that changes with C„, with the poloidal-field 
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Table 3. Ratio of the final intensity of the magnetic field in con- 
figurations A and B. The values in the table are calculated for 
1^1 > 10. For 1^1 < 10 we still have Sf"(A)/Sf"(B) > 1 and 
S|j"(A)/Sp"(B) < 1, but the exact value depends on the differen- 
tial rotation. 
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13.0 
0.37 



ratio being essentially independent of the differential rotation or 
the spin period. 

4.4. Asymptotic states of the static configurations 

The value of the magnetic field after 40 s of evolution for some 
representative cases of the static configurations is shown in Figs. 
|6]-|9] which report the final intensity of the toroidal and poloidal 
components as functions of \q\ and for different values of C„. 

Figures |6] and |2] refer to configurations without //-quenching 
and show that the qualitative evolution of the field is independent 
of the rotation rate of the star. Furthermore, for small values of 
the differential rotation, i.e., for \q\ < \q*\, where |^*| is a repre- 
sentative threshold value, both the toroidal S, and poloidal Sp 
components of the magnetic field are constant with \q\, while for 
1^1 > 1^*1, S, begins to increase and Sp to decrease, following 
power-laws with exponents of the same magnitude but opposite 
signs (see Table|4]i. Note that the value of \q*\ is smaller for high 
spin rates since the dynamo is more easily excited when the star 
is rapidly rotating. 
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Fig. 6. Intensity of the final toroidal field S, in units of Bgq, at the 
end of the 40 s NFI evolution time, as a function of the differen- 
tial rotation parameter \q\ for different values of €„■ The simula- 
tion refers to the configuration A, thus without 77-quenching. 
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Fig. 7. The same as in Fig. |6]but for Sp 



Figures [8]-|9] on the other hand, show the influence of 77- 
quenching on the mean-field dynamo by comparing the final in- 
tensities of the magnetic fields for configurations for which the 
quenching is either active or not, i.e., configurations Aq and A, 
respectively. It is quite evident that the main consequence of rj- 
quenching is to increase the amplification factor of the dynamo 
by several orders of magnitude, leaving the qualitative behaviour 
of the magnetic field unchanged, except for a value of \q*\ sig- 
nificantly smaller. 

A way of interpreting these results is to recall that the inclu- 
sion of //-quenching enhances the transformation of the poloidal 
magnetic field into a toroidal one (i.e., fi-effect) and therefore 
enhances the importance of differential rotation. In addition, by 
reducing the magnetic diffusion, //-quenching effectively favors 
the amplification of the magnetic field (which becomes essen- 
tially frozen with the fluid), increasing the exponential growth- 
rate of the intensity; in the case shown in Fig. [TT]it goes from 



' amp 



0.6 td to T, 



amp 



0.8 Td. After a few diffusion times the 
a-quenching stops the amplification, leading the system towards 
saturation. 

A typical example of the spatial distribution behaviour of 
the magnetic field is reported in Fig.[TOl where we show the be- 
haviour of the toroidal field strength, in terms of Beq, as a func- 
tion of the z coordinate at different time steps for the configu- 
ration Aq. As expected, the field is mainly localized in the NFI 
regions where a ^ 0. 

A couple of remarks are worth making at this point. Firstly, 
while the figures in this section refer only to the configurations 
A and Aq, the same behaviour has also been found for the con- 
figurations B and Bq, thus indicating that the results presented 
here are independent of the extent of the NFI zone. Secondly, 
the extreme magnetic-field amplifications obtained are a direct 
consequence of the idealized setup used in our modeling of the 
PNS. We expect that a more realistic description of the geome- 
try of the star and a consistent treatment of the feedback of the 
magnetic field on the dynamics of the plasma will lower these 
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Fig. 8. Intensity of the toroidal field S, in units of Bgq, at the end 
of the 40 s NFI evolution time, plotted as a function of the differ- 
ential rotation parameter \q\ for C„ - 6 (167 ms < P < 500 ms), 
and an initial seed magnetic field of 10"^ Beq. The simulation 
shows the comparison between a configuration with 77-quenching 
(Aq) and without //-quenching (A). 



Fig. 10. Spatial distribution of the toroidal component of the 
magnetic field at different time steps, for configuration Aq, with 
Ca = 10 and \q\ = 0.1. For this configuration a is different from 
zero only in regions with |z| > 0.3. 
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Fig. 9. The same as in the Fig. |8]but for Sp. 



Fig. 11. Behaviour of the 2-norm of Sp as a function of the num- 
ber of diffusion times (f/r) for the A and Aq configurations, 
\q\ - 2, Ca - 4 (250 ms < P < 750 ms), and for seed mag- 
netic field strengths of Bs = 10 "^Beq- The growth-rate of the 
magnetic field during its first amplification phase is the same for 
both configurations, which are eventually driven to saturation by 
the a-quenching . 



estimates, yielding magnetic fields which are less strong but still 
comparable or larger than the equipartition value. 
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4.5. Asymptotic states of the dynamical configurations 

Although equally idealized, the two dynamical configurations 
AB and ABq are expected to provide a better modeling of the 
first 40 s of life of the PNS, during which the thickness of the 
NFI region is assumed to vary from about 70% of the PNS ra- 
dius, down to about to 20% of it (cf. the profiles of a and rj in 
Fig.|2|. However, besides the fact that for these configurations it 
is not possible to define a unique critical period, since it depends 
on the NFI zone thickness {cf. Sect. 14. 2t . we have found that the 
qualitative behaviour of the magnetic field for dynamical config- 
urations is very similar to that discussed in the previous section 
for the static configurations. 

More specifically the most salient difference is that, because 
of the shrinking of the instability zone, the overall amplification 
is reduced with respect to the A configurations and increased 
with respect to the B configurations. Indeed, at the beginning of 
the evolution the NFI zone is as large as that of the A configu- 
rations, and at the end it is as large as that of the B configura- 
tions. Nevertheless, if one considers a spin rate and a differential 
rotation strength high enough to have the dynamo mechanism 
active during all the instability period {i.e. €„ and \q\ larger than 
the critical values of the B configurations), then the final field is 
comparable to or even larger than, that in the static A configura- 
tions. 



"Elsasser" number A = B^I%kqQ.v (Riidiger & HoUerbach 
I2004I I, namely the ratio of magnetic to zonal flow energy. If 
A «; I the backreaction effect is negligible, while for A s> 1 
it is of overwhelming importance, and should produce a strong 
reduction of differential rotation. If we define an equipartition 
Aeq - BI l%nQQ.ri() and assume that the turbulent viscosity v is 
of the same order of turbulent magnetic diffusion (i.e., both mo- 
mentum and magnetic field are transported by the same eddies), 
we have A = /i^Agq, where h - B,IB^a,- For the typical values 
of our PNSs given in Sect. |2] A ^ 0.05 at equipartition, while 
for magnetic field strengths exceeding the equipartition value by 
two or three orders of magnitude, the effect of the Lorentz force 
on differential rotation cannot be neglected. 

Our simplified model does not allow us to predict the inten- 
sity of this effect, even though we obviously expect that the final 
magnetic field will be comparable to or slightly larger than the 
equipartion one. It is also worth pointing out that a reduction or 
suppression of differential rotation does not necessarily imply a 
suppression of the dynamo, since the a^Q can shift to a pure 
a^ dynamo. Furthermore we cannot either predict whether the 
strong toroidal fields generated are stable against the magneto- 
rotational or Tayler instabilities, which would require a dynami- 
cal analysis, beyond the scope of the present exploratory study; 
this will be addressed in a more realistic model we intend to de- 
velop in the future. 



4.6. Tine effect of Lorentz force backreaction on differentiai 
rotation 



Several authors (e. g., We iss et al. 119841 Belvedere et al. 119901 
Roald & Thomas 119971 Moss & Brooke [2003 Oilman & 
Rempel 2005] Covas et al. 120051 Rempel 120061) investigated 
the effect of the backreaction generated by the Lorentz-force on 
plasma motion; this is sometimes referred to as the "Malkus- 
Proctor effect". These investigations indicate that the dynamo- 
intensified toroidal field interacts, via the Lorentz force, with the 
zonal flow (predominantly azimuthal), thus limiting the growth 
of differential rotation and reducing its strength. An obvious con- 
sequence of this effect is that the amplification of the toroidal 
field itself is also diminished. All of the above mentioned stud- 
ies, which were mainly devoted to the analysis of the Sun or 
solar-type stars, made use of an additional Euler equation in or- 
der to include the Lorentz force consistently. 

This is different from what is done in our simple, idealized 
model case, where we use the induction equation only, and as- 
sume the rotational stress AO to remain constant in time for each 
configuration. As a result, any phenomenological, parameterized 
effect of the magnetic field on the akeady fixed AQ would not 
be consistent with our basic assumptions. 

Nevertheless we performed some runs with a phenomeno- 
logical quenching function applied to Cn (Q-quenching), of the 
same type as those applied to a and rj. The result of these 
runs shows a moderate reduction of the final strength of the 
toroidal field for all the examined cases. This result, however, 
might be somewhat misleading. Our dynamo is essentially an 
a^Q. one, which can shift to an aO. dynamo at strong rota- 
tional stresses, with dynamo numbers (see Sect. 14.7b given by 
the products C^Cn and C„Cn respectively. Therefore, the use 
of an Q-quenching applied to Cn essentially amounts to apply- 
ing a stronger quenching to a (and this explains the reduction of 
the toroidal field strength), but does not limit the fixed rotational 
stress. 

A rough estimate of the importance of magnetic field back- 
reaction on differential rotation can be made in terms of the 



4. 7. A generai expression for fhe finai magnetic fieid 

From the discussion about the time evolution of the solutions 
of Eqs. ( [Tol l and (fTTT i {cf. Sect. 14.2! ) and about their asymptotic 
states {cf. Sect. l4.4l and l4.5l l. it is apparent that the configurations 
considered in this work have common features that we believe 
reflect a fundamental behaviour of the mean field dynamo pro- 
cess. Thus we expect these features to be present also when the 
toy model considered here is replaced by a more realistic one. In 
what follows we discuss how to summarize these analogies by 
presenting a general expression for the final magnetic field. 

Using Figs. |6] and |7] as guides, it is easy to recognize the 
existence of a transition value \q*\ such that for \q\ < \q*\ the 
final magnetic field does not depend on the degree of differen- 
tial rotation, while for \q\ > \q*\ it increases as a power-law. We 
therefore express this increase with a phenomenological relation 
of the type 






(21) 
(22) 



where the indices f and p refer to the toroidal and poloidal com- 
ponents, respectively. Note that the threshold value |^*| may well 
differ between ( |2TI ) and (|22]) . and that the constants K,p depend 
on the particular configuration, i.e. on the profile of a and 77 and 
on 77-quenching. 

Using a power-law fit, we have calculated the values for the 
exponents in the phenomenological expressions (1211 1 and ( |22] ) 
and collected them in Table |4] for both the static and the dynam- 
ical models, with and without 77-quenching. The reported values 
of ji^p have been computed using several configurations having 
different Ca (between 6 and 200), while those for 6,^p have been 
derived from configurations differing by the amount of differen- 
tial rotation {\q\ in the range 10"' - 10^). For all of the configura- 
tions considered, the variance around the reported values is very 
small and of a few percent only. 

Overall, the data in Table |4] show that both -y, and jp are 
very close to either 1/2 or -1/2 (with the exception of Aq, for 
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Table 4. Indices of the power-law behaviour found in the depen- 
dence of B*'" on the degree of differential rotation \q\ and the spin 
rate €„■ The errors, not reported in the table, are of a few percent 
only. 
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which they are ~ 0.75 and ~ -0.35, respectively) and that quite 
generically |y,| + \yp\ - 1. We can therefore rewrite expressions 
(l2Tl i- (l22] | simply as 



^fin 



P 



Kp \q\-"^ . 



(23) 
(24) 



The behaviour of the toroidal and poloidal final fields in the 
super-critical dynamo regime, as given by Eqs. ( l23T l and (l24l i. can 
be understood in terms of the characteristic dynamo parameters, 
Cq, and Cq. The first parameter represents the dynamo's ability to 
regenerate the poloidal field from the toroidal one (i.e., or-effect), 
while the second parameter the dynamo's ability to regenerate 
the toroidal magnetic field from its parent poloidal field, via the 
action of differential rotation (i.e., Q-effect). Both parameters are 
important as they define the "dynamo number", given by Naa = 
CaCn for an aQ. dynamo, and N^^q - C,jCn for an a^Q. dynamo. 
In essence, below a certain critical number A^c, which depends 
on the specific model considered, no dynamo action is possible, 
while the dynamo becomes ever more efficient as A^ increases. 

In our case, when the effect of differential rotation becomes 
important (i.e., high \q\ values), our dynamo behaves essentially 
as an aQ. dynamo with Cn/Ca - 5 |^| » 1 [see Eq. (fT3bl. 
For this dynamo, the ratio of the toroidal to the poloidal field, 
St/Sp - {CnlCaY^^, is high, as verified also with our nonlin- 
ear calculations (see Figs.|6]and|7]l, with typical values ranging 
from 20 to 30. The relationships described by Eqs. ( |23] ) and ( l24l l 
have been derived in the super-critical regime region where S, 
substantially increases and Sp decreases with increasing Cn or 
\q\ for each fixed Ca value (see Figs. |6] and |7] for log|^| > 0). 
Therefore, it is not surprising that S, oc A'^^ [Eq. (l23Tl1, namely 
that the toroidal field increases with dynamo efficiency. On the 
other hand, the poloidal field Sp oc CJN^^^ [Eq. (|24|] will nec- 
essarily decrease with dynamo efficiency increase, since the in- 
crease of A is only due to the increase of Ca (C„ is kept fixed 
for each simulation). This is consistent with the fact that, at 
high differential rotation rates, the zonal flow dominates with 
respect to the vertical motions (predominantly radial) thus re- 
ducing the strength of convection and therefore the efficiency of 
the (T-effect, namely the regeneration of the poloidal field. 

We expect that more sophisticated calculations will produce 
changes of these values reported in Eqs. (|23] | and (l24l i. espe- 
cially in the transition between the a^Q and the orfi dynamos, 
and when the Lorentz-force feedback is properly taken into ac- 
count. However, we expect that the behaviour of the magnetic 
field as given by Eqs. ( |23] l and ( l24l l will remain unaltered. 



5. Conclusions 

We have presented a toy model to describe the amplification of 
the magnetic field inside a proto-neutron star (PNS) via a dy- 
namo action. The model assumes that a neutron-finger instabil- 
ity (NFI) develops in the outer regions of a PNS during the early 
stages of its life as discussed by Miralles et al. ( 200011200? ). and 
that the conditions for the generation of a mean-field dynamo 
process are met. Although highly simplified by being only one- 
dimensional and by adopting the kinematic approximation, our 
model aims to capture the qualitative features of the dynamo ac- 
tion by including a moving boundary of the instability zone and 
the nonlinearities introduced by the feedback processes, which 
saturate the growth of the magnetic field (i.e. a-quenching) and 
suppress its turbulent diffusion (i.e. //-quenching). 

In essence, the amplification of the magnetic field is de- 
scribed in terms of a system of coupled partial differential equa- 
tions of mixed hyperbolic-parabolic type, which are solved nu- 
merically for a very large variety of initial conditions. These in- 
clude varying the spin period of the PNS, the strength of the 
differential rotation between the core and the surface, the inten- 
sity of the primordial (seed) magnetic field, and the extent of the 
NFI zone. 

Overall, we have found that, independently of whether the 
size of the NFI zone varies in time or not, the amplification 
of the magnetic field undergoes a first exponential increase 
with growth-time that is the same for both the toroidal and the 
poloidal components of the magnetic field. The exact value of 
the growth-time depends on several parameters and it is roughly 
in the range Tgmp ~ [0.5 - 2.5] tq, with To being the diffusion 
timescale. The exponential growth then stops through the back 
reaction of the c-quenching and the magnetic field reaches sat- 
uration. The final magnetic field produced at the end of the 40 s 
of evolution does not depend sensitively on the initial magnetic 
field, but it does depend on whether the //-quenching is active or 
not, becoming 2-3 orders of magnitude larger in the first case. 

Despite its crudeness, our model is also able to capture an- 
other important feature of ffie dynamo mechanism, namely the 
existence of a critical rotation period P^, above which no dy- 
namo action is possible and the magnetic field simply decays 
(Bonanno et al. ,2003). For periods near ffie critical one, on 
the other hand, the dynamo is just able to sustain the magnetic 
field close to its initial value, thus avoiding its decay. However, 
as the spin rate (or the degree of differential rotation) is in- 
creased, the dynamo becomes more and more efficient, amplify- 
ing the magnetic field up to values several orders of magnitude 
larger than the equipartion magnetic field. These very high in- 
tensities (10'^ G for the toroidal component and 10'^ G for the 
poloidal one) may seem unphysical at first sight; however, here 
we are considering magnetic fields still inside the neutron star 
and only the poloidal component is thought to emerge after- 
wards. Determining the critical period accurately is important 
to constrain the fraction of neutron stars that may undergo this 
magnetic-field amplification at birth, and we have found that Pc 
is in the range 33 - 600 ms for rigidly rotating PNSs, becoming 
larger as the degree of differential rotation is increased. As a re- 
sult, as long as we have no general constraints on the strength of 
the differential rotation, a lower limit of P ^ 30 ms can be taken 
as the generic threshold below which a mean-field dynamo may 
be active in a newly born PNS. 

Another interesting result of this investigation is that, despite 
the large parameter space considered, the final value of the mag- 
netic field seems to follow a surprisingly robust dependence with 
the spin period and the degree of differential rotation, both of 
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which can be summarized in a phenomenological expression of 
the type 

Sf'" ex {Cj\qV (25) 

which holds only for \q\ larger than a transition level of differen- 
tial rotation \q'\, whose exact value depends on the configuration 
and on the other parameters (overall it is in the range 10~^-10'). 
The exponents 6, y are different for the toroidal and poloidal 
magnetic field components and depend only very weakly on all 
of the parameters varied in this analysis. In particular, for the 
toroidal magnetic-field component we have found 6 ~ \ and 
y ~ 1/2, while for the poloidal one 6 ~ Q and y 1/2. The ex- 
act values of these exponents are likely to be modified by more 
realistic and multidimensional calculations, but we also expect 
that the scaling in expression ( |25] | will persist in further refine- 
ments of this treatment. 

The work presented here can be improved in a number of 
different ways. A first possibility is a more realistic description 
of the geometry of the problem, with a two or three-dimensional 
description of the PNS. A second and computationally less ex- 
pensive alternative is that of improving the nonlinear feedback 
of the magnetic field on the dynamics of the matter This can be 
done by using the same geometry adopted here, but coupling the 
mean-field induction equation ([T]i with the solution of the MHD 
equations for the conservation of energy and momentum. 
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